Global stability for SIR and SIRS models with differential mortality 
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Abstract 



We consider SIR and SIRS models with difTerential mortality. Global stability of equilibria is established 
by using Lyapunov's method. 
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1 Introduction 

The SIR model is a classical model in mathematical epidemiology. Particularly Kermack and McKendrick 
[5n| use a SIR model to prove the existence of threshold. The model of Kermack and McKendrick is without 
demography, i.e. without vital dynamics. The classic SIR models are very important as conceptual models 
(similar to predator-prey and competing species models in ecology). The SIR epidemic modeling yields the 
useful concept of the threshold. 

When a vital dynamic is introduced the asymptotic behavior changes. When the death rates are equal 
in each compartments S,I and R, and equal to the birth rate, the global stability has been solved in 
Actually, the total population is constant, hence the system reduces to a two dimensional system. Then 
using phase plane methods (Poincare-Birkhoff) and Lyapunov functions the global stability is obtained. 

Models with a variable total population size are often more difficult to analyze mathematically because 
the population size is an additional variable which is governed by a differential equation 

The global stability using Lyapunov functions of SIR model with a total constant population is proved 
in [10 ) . However in this model the death rates of S and / are equal and the death rate of the removed 
compartment is adjusted relatively to the death rate of S and the constant birth rate. This adjustment 
is just done to have a constant total population. This is a little bit artificial. The model with constant 
population simplifies in two important ways : 

SI 

• The mass action law — reduces to a bilinear law j3S I 

• The system is a two-dimensional system. 

In this paper we propose a more realistic model, with constant population. We suppose, which is more 
or less observable, that the natality compensates for the mortality. Our model can deal with different death 
rates, and particularly with a over-mortality from the disease. 

We denote TZq the basic reproduction number. It is defined as the expected number of new cases of 
infection caused by a typical infected in a population susceptible [3~ll 16j . We prove in this paper the global 
stability of disease free equilibrium (DFE) if TZo < 1 and that there exists a unique endemic equilibrium 
(EE) if TZo > 1, which is globally asymptotically stable on the domain minus the stable manifold of the 



The stability analysis of classical SIR model is well know since 1976 [53 |T]. The reason was that study of 
stability for these models reduce to the study of 2-dimensionnal systems, hence phase methods can be used: 
Poincare-Bendixon theorem. Periodic orbits are ruled out using Dulac criteria or a condition of Busenberg 
and Van Den Driessche ["2] . 

In the recent litterature, the Lyapunov method is successfuly used to prove the global stability of endemic 
equilibrium. This method consists to find one function, called Lyapunov function and usually denoted by V, 
positive definite and its derivative along trajectories is negative definite. If the derivative V is only negative , 
the LaSalle's invariance principle extend the Lyapunov method in particular cases. This Lyapunov function 
is very difficult to exhibit. However, the class of Lyapunov function 
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is used. This function has a long history of apphcation to Lotka-Voltera models and was discovered by 
Voltera himself, although he did not use the vocabulary and the theory of Lyapunov functions. In 2002, 
Korobeinikov and Wake use this type of function to prove the global stability for SIR, SIRS and 575 
models [10 | and in 2004, for SEIR and SEIS model [5] and give a simply proof of the result of Li and 
Muldowney [15 ) . 

We give a brief outline of the paper. In section 2, we formulate the model and study the global stability 
of the DFE if TZo < 1. In section 3, we study the global stability of endemic equilibrium if TZq > 1. In 
section 4, we deduce results of the next section to SIRS model. In finally we conclude in section 5. 



2 Model formulation 

We consider a population A*' divided into classes of susceptible, infectious and removed individuals, with 
numbers at time t denoted by S(t), I(t) and R(t) respectively, that is A'^ = S{t) + I{t) + R{t). We assume 
that there is no vertical transmission, then all offsprings are susceptibles. We assume that the natality A 
compensates for the deaths. Then A = niS + fi2l + l-iaR- The parameter 7 is the rate of recovery. Note 
that, in our model the disease confers a permanent immunity. The parameter j3 is the effective per capita 
contact rate of infective individuals. We modelize the contact by the classical law of mass action mass. We 
have the following flow graph: 




The dynamic of this model is given by the following system: 

r 5 = A-/?^-^ii5 

{ / = /?^_^2/-^J (1) 

[ R^jI-fisR 

Which reduces to 

5 = -p^ + ^l2I + ^i■^R 

j = /?^-M2/-7/ (2) 
7? = 7/-At3i? 

The population size is constant, so that 5+/ + -R = A'', then we can omit the equation of removed population. 
We obtain the two-dimensional system: 

r S = -l3^+n2l + ti-A{N-S--I) 

For simplicity we can consider the prevalence, i.e. the proportions. 

If we denote j^, the susceptible and infectious fractions, again by 5 and /. Then the system ([2]) is 
reduced to 

J S ^ fi-s + {fJ.2 - tJ--s)I - H'sS - PSI , . 

1 i = psI-{^^2 + J)I ^ ' 

We have < 5, < / and S + 1 < 1. The biological domain of this two-dimensional system is the standard 
simplex. 
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The set Q, = {{S, 1) : S > 0; / > 0; 5 + / < 1} is a positively invariant compact set for The system is 
well posed. 

The basic reproduction ratio is given by 



7^o = 



+ 7 



2.1 stability of DFE 

System ((3]) has a disease free equilibrium state, which is given by (5'*,0) = (1,0) . 

Theorem 2.1 If TZo < 1 then the DFE is globally asymptotically stable on Q,. 
Proof: 

We consider the Lyapunov-LaSalle function V{S,I) = /. We have: 

V = i 

= PSI - (m2 + 7) ^ 
= 7(7eoS-l)(M2+7) 
< 

Furthermore V = if J = or 5 = S* and TZo ~ 1- Hence the largest invariant set contained in the set 

is reduced to the DFE. Since we are in a compact positively invariant set, 
by the LaSalle's Invariance Principle |12 ifT] , the DFE is globaly asymtotically stable in 57. 

Remark 2.1 Unlike Lyapunov's theorems, LaSalle's principle does not require the function V{x) to be 
positive definite . If the largest invariant set M , contained in the set E of points where V vanishes, is reduced 
to the equilibrium point, i.e. if M — {xo}, the LaSalle's principle allows to conlude that the equilibrium is 
attractive. But a drawback of Lasalle 's principle, when significant, is that it proves only the attractivity of the 
equilibrium point. It is well known that in the nonlinear case attractivity does not imply stability. But when 
the function V is not positive definite, Lyapunov stability must be proven. This is why LaSalle's principle 
is often misquoted. Some additional condition enables, with LaSalle's principle, to ascertain asymptotic 
stability. To obtain stability from LaSalle's principle some additional work is needed. The most complete 
results, in the direction of Lasalle's principle to prove asymptotic stability, have been obtained by LaSalle 
himself (haSa\\e:ll3 j. m 1968, completed m 1976 |I7J 

3 Global Stability of endemic equilibrium 

An equilibrium for system ((Sj, different from the DFE, is given by (S, /), where 

S=tl+2^1- and I^^^(l-± 
P 7^o MS + 7 V •^0 

This equilibrium is in the simplex, i.e., < 5", < / and S + / < 1 iff TZo > 1- 
Clearly < 7 is equivalent to 72,o > 1. Now we can write 



1^ +/^3 



S + I = 



7 + M3 

When TZo = 1 this equilibrium coincides with the DFE. Then there is an unique equilibrium in the interior 
of the simplex iff TZo > 1. 

Theorem 3.1 If TZo > 1, the DFE is unstable and there exists a unique endemic equilibrium {S,I) and this 
endemic equilibrium is globally asymptotically stable on the domain Q \ [0, 1] x {0}. In other words on the 
simplex minus the stable manifold of the DFE 

Proof: 

When TZo > 1 the instability of the DFE comes from 3 . 

Let Qi the set defined by = 1(5",/)/ 5" > ~ , 7 > 0, S' + 7 < l|. The set Qi is a compact 

O 

positively invariant. We Consider on fli the Lyapunov function defined by 

v{s, i) = {s-s)^t^ log -^^ + ^- + ^^ + (/-/)- / log i 

P -M2 + M3 + 7 
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It is easy to verifiy that V is definite positive, that is V{S, I) > and V{S, /) = if and only if {S, I) = {S, I). 
His derivative along trajectories of Q is given by: 



V{S,I) 



b - (IJ.3+ J) ■ h 



/i2 + MS + I^S 

/3S/-(m2+7)^-/(/35-(m2 + 7)) 



+ tJ,3 + I3S 



/?S/-(m2+7)/-/(/?S-(m2 + 7)) 



(1 - 5) - (m3 + 7) ..^':^L - -f (/SS - (m2 + 7)) 



M3 (1 - 5-) 
M3 (1 - 5") 

(S - 5) 



-fJ-2 + IJ.3 + I3S 
fJ-3 +'7 



+ M3 + /3S 
-/3S + /JS* 



/(/3S'-(M2+7)) 



Ai2 + M3 + PS 

1 - s 



M3 



M3 +7 
1 - 5" 



(1 - 5) - I3S) 



— At2 + M3 + /SS A13+7. 
1-S 1-5 



-IJ.2 + fJ.3 + PS -^i2+M3+/35' 



/?M3 

At3 + 7 



-M2 + /3 + M3 



-At2 + 



< 



Then we conclude V is semi-definite positive. Then the endemic equilibrium is stable by Lyapunov theorems. 
We prove the attractivity of endemic equilibrium using Lasalle's principle. 

The set on which V — is given by i5 = 1(5*, /) Gf2i / S = s| . Then on this set, we have 5* = /13 + (/12 — 
fj,^)! — IJ.3S — PSI = 0, then / = ^^^^^ — = /. Furthermore the largest invariant set contained in the 



set 



PS - H2+ M3 

I (5*, /) EQi I V {S, I) = o| is reduced to the endemic equilibrium. Hence {S, I) is attractive. Then EE 



is GAS on Qi. 

If 5 < ^^2^-^, we have: 

5 = m + {fj.2 — fJ-a)! - /J-aS — I3SI 
= M3 (1 - 5-) + (m2 - M3 - / 
> 

00 o o 

Then S > 0. Furthermore all trajectories in SI \ fii enter in Sli. Then the set Q,i is absorbant. Hence the 
EE is GAS on Q, . 

In the boundary S = 0et5 + /=1, the vector field is strictly pointing inside Q. Only the S-axis is 
invariant. The endemic equilibrium is GAS on f2\ {{S, I) : I = 0;0 < S < 1}. This end the proof. 



4 SIRS Model 

In this section, we consider a SIRS model with different mortality. With the same notation that above, we 
have the following system: 



Which reduces to 



S^^P^+fl2l+{^i3+l^)R 

^ = /3¥-(M2+7)/ 



(4) 



(5) 



I, R^jl - {ti3 + '^)R 
The system ([5]) is exactly as system ((2)1 where jj,3 is replaced hy (13 + v. 

5 Conclusion 

In this contribution, we have proved the global stability of SIR and SIRS models with differential mortality 
by Lyapunov methods. Our results encompass and improve the results of [10 ] . 
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